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Heat Transfer Analysis of a Closed Brayton Cycle Space Radiator 

Albert J. Juhasz 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

This paper presents a mathematical analysis of the heat transfer processes taking place in a radiator 
for a closed cycle gas turbine (CCGT), also referred to as a Closed Brayton Cycle (CBC) space power 
system. The resulting equations and relationships have been incorporated into a radiator sub-routine of a 
numerical triple objective CCGT optimization program to determine operating conditions yielding 
maximum cycle efficiency, minimum radiator area and minimum overall systems mass. Study results 
should be of interest to numerical modeling of closed cycle Brayton space power systems and to the 
design of fluid cooled radiators in general. 


Introduction 

Due to their relatively low mean effective heat rejection temperature resulting from their 
characteristic “near isobaric” heat rejection process, Closed Brayton Cycle (CBC) space power systems 
tend to require large radiator surface areas for a specified heat rejection quantity, as dictated by the cycle 
temperature ratio and working fluid mass flow rate and the specified by power level. Consequently the 
radiator and associated heat rejection subsystem accounts for a major fraction of the overall CBC power 
system mass. Hence a mathematical analysis of the multi- variable heat transfer processes which affect 
radiator area requirements should provide an important analytical tool in minimizing radiator area and 
thus arrive at an optimized CBC configuration for given cycle operating conditions. 

It is for the above mentioned reasons that a detailed one dimensional analysis of the interacting 
conduction- convection-radiation heat transfer processes was undertaken. The analysis focuses on the 
CBC working fluid heat rejection components, consisting of a fluid transport duct penetrated by 
evaporator sections of heat pipes (HP) which absorb heat from a heat transport fluid flowing over them by 
convection heat transfer. The finned condenser parts of these HP are located outside the fluid transport 
duct and are thus exposed to the space environment. It is these HP condenser sections that transfer the 
cycle reject heat by conduction to high thermal conductivity fins which accomplish the final step in the 
heat rejection process by thermal radiation to the space environment whose equilibrium sink temperature 
can be determined (Juhasz 2001). 

A literature search of CBC space power systems revealed only sporadic reference to CBC radiator 
area computation under certain simplifying assumptions. Klann (1970) for example shows an equation for 
computing radiator area under the assumption of infinite internal heat transfer. This assumption sets 
internal wall surface temperatures equal to the higher gas temperatures. Hence the radiator average outer 
wall temperatures are also computed to be higher than their true values resulting in radiator areas that are 
too low. Hence the analysis of this paper was undertaken to document the derivation of pertinent 
equations, using conservation of mass and energy. Note that since the analysis is limited to determining 
radiator area requirements for specified mass flow rates and temperatures, but not radiator pressure drop, 
only conservation of mass and energy relationships were considered. Since the momentum equation was 
not required for the determination of radiator area, it is not included in the analysis. The results are in 
agreement with earlier work by the author (Juhasz 2005) which integrates radiator heat rejection with 
CBC performance and mass analysis. 
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Applicability of Analysis to CBC Configurations 

The mathematical derivations for analyzing the space heat rejection process, including the 
determination of total radiator area and wall surface temperature profile are valid for various CBC 
configurations as suggested by figures 1 and 2. The relationship of the direct heat rejection radiator for a 
non-regenerated system with direct heating by means of a high temperature gas reactor is illustrated by 
the cycle schematic of figure 1. Cycle heat rejection occurs by the gaseous working fluid flowing over 
evaporator sections of heat pipes whose finned condenser sections radiate to space. In contrast the indirect 
heat rejection scheme shown in figure 2 comprises transfer of heat by means of a heat exchanger from the 
gaseous cycle working fluid to an arbitrary radiator coolant fluid which is a pumped liquid. 

The heat transfer within the radiator is now from the liquid radiator coolant flowing over the heat pipe 
evaporators and the final heat transfer step from the finned heat pipe condensers by thermal radiation is 
analogous to the “direct heat rejection” case of figure 1. 


REACTOR POWER HEAT 

HEAT SOURCE CONVERSION REJECTION 



Figure 1. — Direct heat rejection radiator for non-regenerated closed Brayton cycle power system. 


REACTOR POWER HEAT 

HEAT SOURCE CONVERSION REJECTION 



Figure 2. — In-direct heat rejection radiator for regenerated closed Brayton cycle power system. 


NASA/TM— 2007-215003 


2 












Derivation of Radiator Energy Balance Equations 

Fluid and Wall Temperatures 

A one-dimensional model of a fluid being cooled while flowing in a closed duct whose outer surface 
rejects heat to the space environment by thermal radiation is suggested by figure 3, which is intended to 
show that a gas, or liquid, cooled by thermal radiation using heat pipes can be analyzed as simple one- 
dimensional diabatic flow in a constant area duct. 

For an elemental radiator length the heat transferred from the fluid to the outer wall must equal the 
heat radiated from the outer wall to the space sink. Hence, we can write: 

hft-TjiAi =u-e - T ;\dA r 


A fluid convective heat transfer coefficient related to the radiating area is defined as 


h r =hj 


dAj 

dA,. 


Using equation (2), equation (1) can be solved for the fluid temperature T, giving: 


r = r +— .(r 4 -t 4 ) 

i w h V w s ) 


Implicit differentiation of (3) results in: 


dT = 


4-ct-s 

v h r 



+ 1 


• dT u 


( 2 ) 


(3) 


(4) 


Twin *■ 


1 wex 



Figure 3. — Analogy of radiation cooled gas flow inside a 1-D duct to 
flow over heat pipe evaporators. 
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Note that with fluid temperature T known eqs. (3) and (4) are used to solve for T w via a Newton- 
Raphson iteration. 

A one-dimensional heat balance over an incremental section of the radiator, dA r expresses the 
equivalence of a liquid or perfect gas enthalpy drop to the thermal radiation heat loss by the same 
incremental area, dA r . 


-m-C p -dT = a-8-(r w 4 -T s 4 ^j-dA r 
Where m is the radiator fluid mass flow rate 


R.y 

C v is the fluid specific heat; where C„ for a perfect gas 

F ' y - 1 

R. 

R is the Gas Constant R = and y is the specific heat ratio R u 

Mol.Wt 

the Universal Gas Constant 


= 8314.34 


Joule 
kmol - K 


is 


Determination of Incremental and Total Radiator Area 

Substituting equation (4) into equation (5) yields the following expression 


- m ■ C , 


f A 'T 3 

4 • CT ■ 8 ■ T ... 


+ 1 


• dT w = a • 8 • ( T w 4 -T s 4 } dA } 


Solving equation (6) for dA r results in 


dA r - -m C P 


4 -r„ 


-KW~ t s 4 ) v*W- t s A \ 


■ dT„ 


The Integral of equation (7) over the entire radiator between the radiator outlet and inlet wall 
temperatures, T wex and T win , respectively, can be expressed as 



Although the first term in brackets of equation (8) can be readily integrated, one needs to resort to the 
method of partial fractions to integrate the second term. The integration steps are 
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which leads to the final expression for A r as shown in equation (10), after substituting the indicated 
integration limits. 


A r —m-Cp ■ 


—Id 


f T 4 _T 4 ^ 
*win 4 5 

T 4 _ T 4 

\*wex 4 s J 


+ 


4- a- e- 71: 


l n ( T win~ T s)\ T wex +T s ) o 
{ rr wex~ T s)\ T win +T s) 


W'-^-tan 4 -^ 


V 


's J 


( 10 ) 


Radiator Area Fraction to Reach a Specified Wall Temperature 

A powerful feature concerning equations (8), (9), and (10) needs to be mentioned. If one is interested 
in determining the fractional radiator area A r (x) to a region where the wall temperature reaches a 

specified value, T w (x) where T we x ^ T w {x)< T win then the equation (9) can be integrated between upper 
and lower limits of T w {x) and Twin 



V (x) 


4 •a-e-T w 

kW~ T s A ) 


1_ 

a-8-(r M 4 



( 11 ) 


Evaluation of the integral shown in equation (11) will result in an expression for /l r ( x j , which is 
similar to equation (10), but for the T wex symbols being replaced by T w (x ) . To determine the fractional 
distance, x, along the radiator where the wall temperature reaches T w (x ) , one can take the ratio of areas to 
obtain: 


In 


f 4 _ 4 N 

win * s 


T w {x) 4 -T s 


4-a-e -TV 


In- 


( T win- T s )-(t w (x)+7;) 

[ T w [ x )~ T s )\ T win +T s ) 


T ■ 

_i x mn 

tan . 


-tan" 


_! TwW" 


V 


(12) 


In 


( 4 ~ 4 h 

*win ^ s 

rji 4 y 4 

ywex ~*S J 


4-a-e-TV 


(■ T win ~ T s )' i T wex + T s) ( 

In 2 

{ T wex~ T s)\ T win +T s) 


T ■ T 

win _i wex 

tan tan 






•5 7 . 


Equation (12) can be used to compute the fraction of total radiator area bounded by specified wall 
temperatures, in order to determine the weight of multi-segment radiators built from different materials as 
dictated by wall temperature requirements. 

A mean effective radiator wall temperature T we jj- can be derived by substituting equation ( 1 0) into a 
heat balance for the entire radiator, namely: 

\ T in - T ex)=VZ- A r\ T M,eff 4 O 3 ) 

where T in and T ex are the working fluid inlet and outlet temperatures, respectively. Solving equation (13) 
for T we a- one obtains: 
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0.25 


(14) 


T weff 


m ■ C i 


( T in ~ T ex ) , T 4 

(<J-S-A f .) S 


All values on the right hand side of equation (14) are known, except for the fluid temperatures T,„ and 
T ex . But these can be expressed in terms of the corresponding wall temperatures by using the relationships 
given in equation (3): 


t =t ■ +— -fr • 4 -r 4 ) (15) 

* in * win 1 ^ y- win s ) V 

and 

T e x=T weX +~{fy V e X A -T s A ) (16) 

Finally substituting equations (A15) and (16) into equation (14) the following expression is obtained 
for T weff 


T weff ~ 

Equation (17) is the “mean effective” — or integrated wall temperature of a duct which is carrying a 
radiatively cooled liquid or gaseous working fluid with known heat transport properties. 

Note that, although the above analysis has been performed with the assumption that the heat transfer 
coefficient, h r , is constant throughout the duct, it is possible to include variable heat transfer conditions as 
discussed next. 


4 m-C 

T 4 + p - 

a-e-A,. 


T ■ -T + — .(r ■ -T ) 

x win ± wex 1 7 V-* win M wex ) 

fl r 


(17) 


The Case of Variable Heat Transfer Coefficients 


For the case where the heat transfer coefficient is a function of temperature and duct position a finite 
difference computational approach was used. The duct was subdivided into n sections with running index 
j ranging from j = 1 to n, with each of these sections having a specified heat transfer coefficient, h,-( j) 
controlling the internal heat transfer process over section j. Equation (11) can then be rewritten as 


r A r ( j + 1 ) 

U-(/) 


dA r = - m C p 


r T w (j+ 1) 
>T w (j) 


4osT~ 


h r (j)\ T w- T s) a-s-(r 4 -r/) 


dT,„ - 


(18) 


The total required heat transfer area may then be found by summing the elemental areas in 
accordance with equation (18) from element j to element n, as shown below. 


A 


tot 


y rMJ+ 1) 
j—' 


dA r 


(19) 
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Equations (18) and (19) were programmed in the radiator subroutine and validated by subdividing the 
radiator into 25 sub-sections each having identical heat transfer coefficients, h r (/ ) . Performing the 
summation indicated by equation (19) the same total area was obtained as for the integral (non-sectioned) 
radiator having the same internal heat transfer coefficient. 

When substituting monotonically increasing heat transfer coefficients in the sectioned radiator 
configuration, the result was a reasonable reduction in total area of ~3 percent. 

Temperature Profile Example 

An illustration of a wall temperature profile resulting from a given radiator fluid temperature 
distribution and with a gradually increasing h as the flow proceeds from the radiator duct inlet to exit will 
be more meaningful if comparative wall temperature profiles with constant h are also shown on the same 
plot as illustrated in figure 4, which shows four “Mathcad” generated curves. The radiator conditions are 
typical for a hypothetical high turbine inlet temperature of 2000 K, with a cycle temperature ratio of 3.5 
and a pressure ratio of 2.0 for a regenerated cycle. 



Fractional Radiator Length 

Figure 4. — Fluid and wall temperature profiles for various internal heat 
transfer coefficients. 

The top curve with the circle symbols, marked Tf, represents a given He gas temperature profile with 
a duct inlet temperature of -860 K and an exit temperature of 567 K. The second curve from the top, T w 9, 
with the diamond symbols shows a corresponding wall temperature profile obtained with a very high 
constant heat transfer coefficient, h = 1300 W/m 2 - K. In contrast, the bottom dashed curve (also with 
diamond symbols, but marked T w ) represents a wall temperature profile for an assumed value of 
h = 100 W/m 1 -K. The third, heavy dashed curve from the top, Twa, depicts the wall temperature 
profile resulting from flow conditions in which the value of h increases monotonically from 1 00 to 
1300 W/m 2 - K. Note that the wall temperature values for this curve agree with the duct entrance value for 
the low and the duct exit value for the high heat transfer case. It is also of interest to note that the variable 
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heat transfer wall temperature profile approaches the wall temperature profile of the high heat transfer 
case asymptotically at the duct exit. This fact provides an important validation check on the correctness of 
the computational procedure inherent in the author’s radiator sub-program. 

A typical example of the influence of variable heat transfer coefficients on the conceptual design of a 
trapezoidal heat pipe radiator for a CBC energy conversion system with a fission reactor heat source is 
discussed in the next section. 


Application of Analysis to Conceptual Design of a Spacecraft Radiator 


In the derivation of the radiator heat balance equations a fundamental assumption was that the 
radiator fluid convective heat transfer coefficient can be related to the radiating area by equation (2), 
namely 


K = K ■ 


dAt_ 

dA r 


( 20 ) 


Concerning the ratio of differential areas, dAj / dA r , it is evident from figure 3 that this ratio is 

identical to the ratio of heat pipe evaporator area inside the duct, which is exposed to the working fluid, to 
the finned condenser area external to the duct, which is exposed to the vacuum of space. 

As an illustrative example figure 5 shows the trapezoidal radiator panels proposed for the SP-100 
conical radiator configuration. Assume that the radiator cooling fluid enters the central duct at the lower 
right of the trapezoidal panel. Furthermore, a constant evaporator length of ~0.2 meters (i.e., for 
100 kWe BC power system) throughout the duct and a gradually decreasing finned HP condenser length 
from 0.6 m at the duct entrance to 0.2 m at the exit are assumed. Thus the incremental area ratio is 
obviously increasing from 0.333 to 1.0. This means that for constant the value of h r will increase by a 
factor of three, even if the duct cross sectional flow area is constant. Of course the opposite trend is 
obtained for radiator working fluid entering the central duct at the short end of the trapezoid panel. In this 
case, for a constant value of the value of h r will decrease by a factor of three. The preferred flow 
direction for minimizing radiator area will be discussed in the next paragraphs using the theoretical 
relationships developed above. 


Heat pipe 
with end caps 


Radiator panel 
with 226 heat pipes 



Design features 
12-panel, conical radiator 
Carbon-carbon heat pipes 
Integral fins 

Potassium working fluid 
Metal liner 


Radiator 
configuration' 


Figure 5. — SP-100 trapezoidal radiator panel and conical radiator assembly. 
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Payload 



Shield 


Reactor 


Radiator 


Turbo 
alternators 


Heat exchangers 
(recuperators) 


Figure 6. — Nuclear-CBC powered spacecraft with trapezoidal HP radiator. 


Of course, for nuclear space power systems the trapezoidal radiator panel shape is dictated by the 
requirement that the radiator fits inside the “shadow” of a radiation shield which is placed between the 
reactor and the spacecraft major components, in order to reduce radiation dose levels (i.e., Neutron and 
gamma) on the power system, radiator, and the payload. Such a conceptual nuclear power system 
comprising a fission reactor with a neutron-gamma “shadow shield,” and a CBC power system with a flat 
panel trapezoidal radiator is shown in figure 6. This figure a conceptual diagram of a 100 kWe class space 
nuclear power system using regenerated CBC conversion technology with the main sub-systems are also 
identified. Note the large trapezoidal area devoted to the radiator which is in the shadow of a conical 
shield of 15° half angle. 

In this system the best upstream point for fluid entry is hypothesized to be at the end where the finned 
heat pipe condenser section is longest, since the heat transfer coefficient, h r would be lowest where the 
radiator fluid temperature is highest. As the radiator fluid proceeds downstream toward the compressor 
inlet, the fluid temperature decreases, but h r is rising, since the ratio dA i / dA r in equation (2) is increasing. 
After the gas has been cooled to compressor inlet temperature at the radiator duct exit, it is piped near 
isothermally in the same direction for re-injection into the compressor. 

To verify the above hypothesis, the derived relationships (especially eq. (18)) were programmed into 
a numerical code which was then used to perform radiator area calculations for the fluid temperature 
profile shown in figure 4, and a total heat load of 1 99 kW t . The computational results for three heat 
transfer coefficient input cases are summarized in table I. 

TABLE I.— RADIATOR PANEL AREA COMPARISONS FOR 3 DIFFERENT INTERNAL HEAT TRANSFER CASES 3 


Case number 

Heat transfer coeff, 
hr, 

W/m 2 - K 

Specific radiator 
panel area, 
nr/kWe 

Radiator area change, 
percent 

Remarks 

1 

200 - Constant 

0.24453 

0 . 

Rectangular 
Radiator panel 

2 

600 to 200 decreasing 

0.22384 

-8.46 

Trapezoidal panel 
Entry @ short end 

3 

200 to 600 increasing 

0.22099 

-9.63 

Trapezoidal panel 
Entry @ long end 


“Total heat rejected =199 kWt; Duct Inlet Temp. = 859.3 K; Duct Exit Temp. = 570 K; Space Sink Temp. = 200 K 


Radiator duct divided into 2 1 segments in flow direction. 
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The following conclusions can be drawn from the three cases shown in this table for the conditions 
indicated by the asterisk: 

(1) For the case of a constant heat transfer coefficient of 200 W/m 2 -K the specific radiator panel area 
was computed at 0.2445 nr/kWe. This value was used as a reference for comparison with cases 2 and 3. 

(2) Case 2 represents a trapezoidal radiator panel with flow entering the radiator duct on the short 
side, as was done for the SP-100 panels. The heat transfer coefficient, h r , was varied from 0.6 to 

0.2 kW/nr K linearly in 21 steps with the direction of the radiator working fluid which in this case was 
identical with the helium gas turbine working fluid. The required panel area is seen to decrease by 
8.46 percent. 

(3) Case 3 represents a trapezoidal radiator panel with flow entering the radiator duct on the long side, 
as is recommended by this author, since this would cause h r to increase from 0.2 to 0.6 kW/m 2 -K. As a 
result the radiator area reduction would further improve from 8.46 to 9.63 percent, as shown. Although 
the improvement over case 2 is only a modest ~1.2 percent, the actual area savings would be significant 
for high power systems requiring several thousand square meters of radiator area. 

The results shown in table 1 show that, as hypothesized, it would be of some advantage to send the 
flow along the central boom to the wide end of the trapezoid, where the finned heat pipe condensers are 
longest and the radiator fluid temperature is highest, and have it enter the duct at this point so that thermal 
energy transfer can occur with gas flowing over the heat pipe evaporators in a direction towards the turbo- 
alternator. 


Conclusions 

A detailed derivation of the mathematical relationships governing heat rejection from a CBC space 
power system was completed. The analysis is valid for both systems having a separate liquid pumped heat 
rejection loop connected to the CBC by means of a gas-to-liquid heat exchanger, or direct heat rejection 
by the cycle working fluid flowing over evaporator sections of heat pipes which protrude into a closed 
heat transport duct. The analysis also showed that designing the radiator duct in a manner that will allow 
the effective heat transfer coefficient to increase in the direction of the flow will result in a modest 
decrease in radiator area requirements. 

Although the analysis was motivated by the heat rejection process of CBC power systems, it can be 
extended to computing surface area and weight of any gas-or liquid cooled space radiator as long as 
required input conditions are known. 


Nomenclature 

Aj inner surface area; A 0 outer surface area 

A f( ■ ) incremental radiator area at section j 

dAj elemental inner wall surface, or heat pipe evaporator area 

dA r dA 0 , is the elemental radiating outer wall surface, or heat pipe condenser area 

R-y 

C v fluid specific heat; where C„ = for a perfect gas 

F F y-1 

F is the view factor to the space environment (F = 1 for this analysis) 

^ heat transfer coefficient 

h 

1 internal heat transfer coefficient 
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heat transfer coefficient normalized to external or radiating area 
radiator fluid mass flow rate 


Gas constant ^ Mol.Wt ) . R u = Universal Gas Constant = 8314.34 Joule/(kmol-K) 

wall temperature at duct outlet 

wall temperature at duct inlet 

is the fluid bulk temperature at axial location x 

is the equilibrium space sink temperature 
T 

u ' (/l wall temperature at arbitrary location x, or j 
effective radiator wall temperature 

fractional distance along radiator flow path 
specific heat ratio of radiator heat transport fluid 

the Stefan-Boltzmann constant = 5.668x10 8 watts/nrK 4 
the emissivity of the radiating surface 
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